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^ . Abstract 

£f~^ ■ We construct a class of simple and calculable theories for the supersymmetry breaking 

soft terms. They are based on quantum modified moduli spaces. These theories do 

not break supersymmetry in their ground state; instead we postulate that we live in a 

■ supersymmetry breaking plateau of false vacua. We demonstrate that tunneling from 

the plateau to the super symmetric ground state is highly supressed. At one loop, the 

Q_i' plateau develops a local minimum which can be anywhere between 10 s GeV and the 

grand unification scale. The value of this minimum is the mass of the messengers of 

supersymmetry breaking. Primordial element abundances indicate that the messengers' 

'"^ ' mass is smaller than 10 12 GeV. 
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1 Introduction 



Gauge mediated theories of supersymmetry breaking [|1| have recently attracted a great deal of 
attention becuse they solve the supersymmetric flavor problem. Beginning with the pioneering 
papers of Dine, Nelson and collaborators 0, a lot of effort has been devoted to the construction 
of explicit realistic models with dynamical supersymmetry breaking @-||. It seems fair to 
say that these models, when they work, are quite involved. 

In this paper we wish to show how to construct simple theories of the supersymmetry 
breaking soft terms. The key ingredient that will allow us to do that is to relax the requirement 
that supersymmetry is spontaneously broken in the full theory. We instead will be content to 
postulate that we live in a supersymmetry breaking local minimum located on a plateau of 
flat directions. 

In section 2 we will present the general class of theories we consider. In section 3 we will 
demonstrate the it is not dangerous to live on a plateau: the tunneling amplitude into the 
supersymmetry preserving true vacuum is infinitesimally small. In section 4 we will present 
explicit models implementing the ideas of section 2. In section 5 we will show that light 
element abundances place a limit that the messenger mass -which is related to the position 
of the local minimum in field space- has to be less than 10 12 GeV. We conclude in section 6 
with some brief comparisons of our models with earlier work. 

2 The General Idea 

We are interested in theories based on a gauge group 

G = SU(5) W x Gb x G s . (1) 

By SU(5)w we mean indeed 577(3) x SU(2) x U{1) as the dynamics we describe takes place 
below the GUT scale. Nonetheless, as we want to preserve perturbative gauge unification, 
we will only consider messenger sectors involving complete SU(5)w multiplets. The role of 
the remaining factors, as it will be clear below, is to provide a stable vacuum with broken 
supersymmetry. We consider the following simple form for the messenger sector superpotential 

W = \ 1 Rqq + \ 2 RQQ (2) 

where R transforms non trivially only under Gb- The q © q are vectorial under SU(5)w, i-e. 
5 © 5, but singlets of Gs- Finally, Q © Q are a vectorial representation of Gs, but singlets 
of SU(5)w All the fields transform non-trivially under Gb- Just to make contact with usual 
gauge mediated models, we anticipate that the q, q will be play the role of the messengers, the 
Q, Q will provide a scale via their condensation, and R will play the role of the usual X field 
getting both a scalar and an F component vacuum expectation value (VEV) from the strong 
dynamics. 
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Furthermore we are interested in theories with the following two properties at the classical 
level: 



I) R contains just one D-flat direction, or equivalently there exist just one holomorphic in- 

variant u = u(R) ~ R k , which is a monomial made out of R elements and k is an integer. 
Along this flat direction the gauge factor G B is broken to a subgroup H B and all elements 
in R but one, which we parametrize by u, are eaten by the SuperHiggs mechanism. For 
notational purposes we will in what follows parametrize the flat direction with a field 
.R, since this is the one with canonical Kahler metric at large field values. 

II) Along u the q, q and Q, Q pair up to get masses which are respectively ~ X±R and ~ A2.R. 



Thus far away along u^O the low energy theory has gauge group G = SU(5)w x H B x G$ and 
the sole matter content is given by the MSSM fields plus the gauge singlet u. The additional 
group factors are just pure gauge with no matter. Furthermore we are going to assume that 
the only sizeable non perturbative low energy dynamics is the one associated with Gs- On a 
case by case basis we will check that H B is either too weak or it does not affect the effective 
superpotential. Above the scale of Gs confinement the effective theory has W — 0. The 
effective holomorphic scale A e ff of the low energy Gs theory is given by the 1-loop matching 
of the gauge couplings in the microscopic and macroscopic theories at the scale X |5) 

A^Js = A^gs-vq^X)^ (3) 

where nc s and hq are the Dynkin indices for the adjoint and Q © Q representations of Gs- 
Below the scale A e // the Gs charges are confined in massive states, and the only effect of 
the Gs dynamics on the low energy theory is to generate a superpotential for X via gaugino 
condensation jlOj [Tl|] 

W eff =< X S X S >= A- 3 ef f = A 3 '^^°s (\ 2 xyQ^Gs . (4) 

In the limit in which the other gauge group factors are treated as weak the above result is 
exact. It can indeed be established by using holomorphy and an i?-symmetry which has no 
Gs anomaly JTIJ. Under this symmetry X has charge 2/i Gs //ig, consistent with the above 



W e ff having charge 2. 

The further requirement that we are going to make in order to generally characterize our 
models is that fiQ = fic s - F° r the class of theories that satisfy this identity W e ff is linear in 
X 

W eff = X 2 A 2 X. (5) 

An example of this possibility is given by Gs = SU(N) with Q + Q corresponding to N 
flavors of fundamentals. This is the crucial property that allows us to build simple and 
realistic theories with locally broken supersymmetry. Eq. (|5]) corresponds in fact to the 
simplest O'Raifertaigh model: F x = dxW e ff = A 2 A 2 , so that supersymmetry is broken at 
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any point on the X complex line. For \lq ^ hg s the scalar potential \dxW e ff\ 2 would push 
X either to the origin or to infinity, where supersymmetry is in general restored. Indeed, in 
models with the simple microscopic superpotential in eq. |2|, at the point X = there are 
in general other flat directions involving the massless q and Q, where supersymmetry can be 
restored. The existence of a plateau V e ff ~ \dxW e ff\ 2 ~ A|A 4 insures that X may stabilize 
in some other way far away from supersymmetry restoring points. This is the novelty of the 
models we discuss: there is no need to break supersymmetry in the true ground state. 

The result in eq. (||) can also be derived by using the superpotential from confinement 
derived in Ref. |13| (see also Ref. fllil ) and by integrating out the massive composite super- 
fields. 

Up to this point we have devised a general class of models for which there exists a direction 
X in field space along which W e ff oc X. In such a situation the curvature of the scalar potential 
is completely determined by the Kahler metric. We are interested in the region X 2 X 3> A, 
where it is consistent to consider only the perturbative corrections to the Kahler potential. 
These we know how to calculate. The tree level Kahler potential is given by (see eg. Ref. 
||15|| ) K tree = (W) 1 ^ = XX^. The radiative corrections are due to loops involving the heavy 
superfields Q,Q and q, q and the heavy vector superfields in the coset Gb/H b . All these 
particles get a mass from the X VEV. At the 1-loop level the result has the form 



K(X, Xt) = XXt 1 + — -{C B g\ - C X X\ - C 2 \ 2 2 ) ln(XXt/M^) + 

107T Z 



(6) 



where the C's are positive Casimirs and the coefficient of the logarithm just corresponds 
to the 1-loop anomalous dimension of the original field R. For large logarithms the higher 
loop terms can be resummed via the Renormalization Group (RG) so that in the leading log 
approximation the Kahler metric reads 

K(X, Xt) = XX^Z R (xxyM 2 Pl ) . (7) 



In the above, Zr represents the R wave function renormalization and X the bare field. Then 
the effective potential is just given by 

Wf A|A^_ 

v "< a x a x ,K ~ z R (xx\) (S> 

where we have approximated dxdx^K with Z(XX^), neglecting small finite corrections. The 
logarithmic evolution of Zr with |X| can generate local minima in the effective potential. 
Indeed the potential energy behaves qualitatively like a squared effective Yukawa coupling. 
As it can be inferred already from eq. (U), in the limit in which g B Ai,A2 we have that 
1/Zr grows at large X thereby preventing a runaway behaviour. Indeed one can see this as 
a signal of the fact that pure Yukawa theories become strong in the ultraviolet by developing 
a Landau pole. On the other hand, in the opposite limit of g B ^> A 1; A 2 the potential grows 
at small X. This can prevent X from reaching the origin or other branches of moduli space 
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where supersymmetry may be restored. It is thus clear that the competition between these 
two effects, which arises when Ai^ ~ Qb m &y stabilize the potential at some local minimum 
Xq. Notice that it is crucial that X comes form a field R which is charged under some gauge 
group Gb- This is precisely what provides the gauge contribution g 2 B without which X is 
always pushed towards the origin, where we loose control of our approximation (and in most 
cases restore supersymmetry). The location of this X is determined by the stationary point 
of the RG equation for Z R 

d\a.\X\ B ^ B ^ D 1 Z R Z q Zq 2 Z r ZqZq 

where gs and the Z's are respectively the running gauge coupling and wave functions (of 
the various fields) evaluated at the scale X, while Ai^ are the bare Yukawa couplings, i.e. 
renormalized at the cut-off scale Mpi. Then the expressions multiplying C\ and C2 represent 
respectively the effective running Yukawa couplings Ai(X) and A 2 (X). This is just the usual 
Coleman- Weinberg mechanism. We have a "perturbative" dimensional transmutation, by 
which it is consistent to expect the minimum X to be in the range A <C X Q <C Mp\. As we 
will explain in detail in the next sections, this is precisely the range which is interesting for 
building realistic models. It is also the range for which the minimum on the plateau, even 
if it is not the true ground state, is nonetheless stable on cosmological time scales (see the 
next section). Finally, to be more precise, in order to insure that the stationary point be 
a minimum one should verify that d 2 Z R / 'rf(lnX) 2 be negative. This equation, see eq. @, 
involves both the gauge /3-function and the 7's of the other fields, and has to be checked case 
by case. 

Before going further we would like to remark that this mechanism for stabilizing a tree level 
flat potential in a supersymmetric theory is just Witten's inverse hierarchy [Tj|. The amusing 
and new thing is that in our case the superpotential is generated non-perturbatively rather 
than being an input at tree level. In other words we have a double dimensional transmutation 

1) < Fx >= A2A 2 is provided by gaugino condensation. 

2 ) < X > is stabilized at some point by the Coleman- Weinberg mechanism. 

It is interesting that in the region X> Awe can determine the exact superpotential together 
with a fairly accurate perturbative approximation to the Kahler potential. This allows us to 
safely control the full scalar lagrangian. We stress once more that, as in Witten's model, it is 
crucial that X be part of a charged multiplet in order to have balancing gauge and Yukawa 
contributions in V e ff. The issue of how a non-perturbative superpotential is stabilized by 
RG evolution is a subtle one, and our result disagrees with ref. ||17|| . The conclusion of that 



paper is that V e ff ~ A 2 A 4 / 'Z r ZqZq, i.e. proportional to the running Yukawa coupling, rather 
than just V e ff ~ A 2 A 4 /Zr which we believe to be the correct result. Notice that the result 



of ref. |17j would allow to stabilize the potential even for a gauge singlet R, as ZqZq already 
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involve a gauge coupling contribution with sign opposite to that of Yukawa's. It seems to us 
that this result violates holomorphy. We believe that the correct approach is to integrate the 
RG without rescaling the chiral superfields. If the fields are rescaled, a spurious dependence 
on A^ and X^ is induced both in the running yukawa coupling and in A e ff, the latter via the 
scale anomaly [f|. In the end, however, this spurious dependence should cancel out and the 
result for W e ff(X) be consistent with holomorphy. We use a scheme in which the fields are 
not rescaled, so that holomorphy is manifest in each step. 

What is the phenomenology of these models? At the decoupling of the messengers q, q, 
effective operators involving the MSSM superfields and X are generated. These are the sources 
of sparticle masses after X gets an F-component VEV. The evaluation of the relevant diagrams 



in component has been performed in ref. [18|. The gaugino masses just correspond to the 1- 



loop contribution of the messengers to the gauge /^-function: 

= J d 2 6N^-ln{\X)W a W a (10) 

where N is the number of messenger flavors and where i = 1,2,3 runs over the 577(3) x 
577(2) x U{\) couplings. The resulting gaugino masses are thus 

M &i F x M ai ( A 2 A 2 \ . . 

m - = N ^ = N ^[z x Jc)- 

where Zx indicates the wave function Zjj(XX^) at the minimum. Phenomenology requires 
AA 2 jZxX ~ 10 — 100 TeV. The sfermion masses arise from 2- loop diagrams, and they are also 
related to the (3 function coefficients that govern the evolution of the chiral superfield wave 
function. For instance the squarks get a contribution oc a 2 and given by 

C Q = J d 4 9N^f (i \n\XXynl v ) + 1h(**7a&v) Mfin/Xtf)} QQ^ (12) 

where \i\jy and fj,m are respectively the ultraviolet and infrared cut-off of the diagram. The 
above expression has a very simple interpretation. The ln 2 (nuv) arises from the integration 
region where both loop momenta are 3> X, while the infrared divergent piece is determined by 
momenta < X in the external loop. This explains the relative 1/2 factor. By expanding the F 
components of X the above expression contributes both to the chiral preserving squark masses 
and to the A-terms. The latter are generated by the simple powers of Fx after eliminating 
the auxilliary quark F-fields. Notice that the contribution to the A-terms is of the form 
a 2 \n(X/ fi IR ), which is precisely the result one would get by solving the 1-loop RG equation 
below the messenger scale. Notice also that it is precisely the relative 1/2 coefficient amont the 
two In 2 contributions in the above equation that ensures that there be no ultraviolet divergent 
contribution to A-terms. This is indeed a constraint one can use to derive eq. (|T2|). Finally 
we report the expression for sfermion masses 

rh 2 = 2 N (^) 2 ±k^r (13) 
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Figure 1: An idealized flat potential with V(0 > A) = U ~ A 4 and V(0) = 0. 

where the sum is over the ST/ (3) x SU(2) x U(l) factors. Moreover k\ = (3/5)(V/2) 2 with 
hypercharge normalized such that Q = T3 + Y/2, k% = 3/4 for SU(2) doublets and zero for 
singlets while k% = 4/3 for color triplets and zero otherwise. 

Since X and Fx are independent scales, we have that their relative ratio can be wherever 
it is consistent with phenomenology and with our approximations. In particular, and we will 
discuss this in more detail later on, the messenger scale X can be anywhere between 10 8 and 
10 14 GeV. Consequently the number of messenger flavors can be bigger than the minimum 
5 which one gets for X ~ 10 5 GeV by requiring perturbative gauge unification |19 
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By varying G$ we can explore essentially all possible gauge mediated boundary conditions on 
sparticle masses. 

We have now to discuss why such a local mimimum on a plateau is cosmologically accept- 
able. 



3 The Stability of our Universe 

The potential V e ff(X) in our class of models can be pictured as a very flat (with logarithmic 
curvature) plateau at X ^> A. Moreover we know that classically at X = there are other 
branches of moduli space. Some of these are supersymmetry preserving even at the quantum 
level. We expect that the transition between the plateau and these other directions of lower 
energy takes place smoothly in the region X < A. Physically what happens in this region, is 
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that some mesons or baryons of Gs become tachionic and develop a VEV. In the Appendix 
we discuss this in some more detail for a specific model. In this section we are going to discuss 
the quantum stability of a false vacuum with these properties. In order to get an estimate of 
the rate of tunneling we can reduce to a one dimensional potential V(4>). In some of the cases 
of interest, the region < A corresponds to more than one vacuum or flat direction. Still our 
conclusions are not qualitatively modified in the realistic case. Furthermore in order to make 
the discussion simpler, we will study the rate of true vacuum bubble nucleation for a truly 
flat plateau, in the limit in which A<X. This is shown in Fig. 1. Indeed, at first sight, this 
may look like a bad approximation as it eliminates the potential barrier that separates the 
physically interesting, but false vacuum, from the supersymmetry preserving one. But this is 
too naive since in field theory the potential energy TC consists of both a gradient and a static 
piece 

U = J c/ 3 x{(V0) 2 + V {<!>)}. (14) 

When the potential V((f>) is flat, it is the gradient energy (V0) 2 that provides the barrier in the 
functional space {</>}. Consider, indeed, a homogeneous field configuration where <f>(x) = X 
over all space. Then for any small and local variation 8(j)(x) of the background field, the static 
energy is unchanged while the gradient contribution is positive definite 

AH = AHgrad + AHstat = AH grad > 0. (15) 

It is also clear that the larger X/A the more gradient energy must be generated in order to 
reach lower potential energy, so the higher the barrier in field space is. We thus expect that for 
large X/ A the rate of quantum tunneling into the true vacuum will be suppressed. The rate of 
bubble nucleation can be then evaluated in the semiclassical approximation. This amounts to 
finding the action for the bounce, a spherically symmetric solution to the Euclidean equations 



of motion p0|. In terms of the bounce action Sb, the rate of bubble formation per unit volume 



T/V is just given by 

£ = ^~ Sb (16) 

where \x is a mass parameter determined by the calculation of quantum fluctuations around 
the bounce and which we expect to be of the order of the largest scale in the problem, X. 
The bounce solution 0(r) satisfies the following boundary conditions and equations of motion 

0(0) = 0(oo) = 0(oo) = X (17) 

+ -0 = ^(0). (18) 
r 

(19) 

The equation of motion describes a classical motion in a one dimensional potential = — V(0), 
with time variable = r and with a friction term whose strength decreases as 1/r. The particle 
analogy is very useful and one can apply considerations similar to those made in Ref. |2(| 



for the thin wall bubble. Similar calculations have been done in Refs. pi, 22\. The size of 
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the bubble is simply established by energy conservation. Consider indeed the motion in the 
plateau region > A, see Fig. 1. Here is just a massless field and the solution to eq. (|T§D 
is simply 

0(r) = (X-A) [l-(^) 21 



+ A. (20) 



Notice that 0(oo) = X and that the radius R is an integration constant characterizing the 
size of the solution: (j)(R) = A. The kinetic energy 2 /2 at r = R is given by (X — A) 2 /2R 2 ~ 
X 2 /2R 2 . In view of the presence of a dissipative force, this has to be smaller than the 
maximum available energy U . Then it must be R > X/y2Uo ^> 1/A, much bigger than the 
typical scale of the potential V(4>). This large radius has two implications: i) as X/A — > oo the 
starting point 0(0) — > 0, where dV = 0, since for 0(0) > 0, where dV ~ A 3 , the radius would 
be too small 0(1/ A), the only scale in the potential; ii) in the region < A the friction term 
are 0(1/ AR) with respect to the leading terms, so that energy conservation holds up to order 
1/RA < A/X. Then from i), ii) and by equating potential and kinetic energy respectively at 
~ and = A we find the bubble radius 

H, = |1<1 + 0<A/*)) *.-|- (21) 

More qualitatively we can describe the bounce solution as follows. starts very close to the 
top, and there it stays until friction is small enough to allow the motion to X ^> A. (Remember 
that in the particle analogy the potential is — V(<p) of Fig. 1.) At r ~ R = X/ \?2Uq the particle 
goes down the hill. This motion is determined by the fundamental scale A in V and happens 
within a thin shell (short time) Ar ~ 1/A <C R- Its contribution to the action is consequently 
negligible. Indeed we can write the action S(R) as 

S(R) = -\tx 2 U^ + 2ir 2 X 2 R 2 + 0(A 3 R 3 ). (22) 
The first term is from the bulk of the bubble < r < R where d> ~ 0. The second comes from 



the massless region r > R, where > A, and by using eq. (20). Finally, the third term is an 



estimate of the contribution from the transition region, the thin shell. By extremizing with 



respect to R we find the same bounce radius as eq. (j21j) . Neglecting the subleading term we 
then get the bounce action 

s *= 2 *1r 2 * 2 © 2 ' (23) 



The consistency of this approximation can also be verified by noticing that eq. Q2lq) is 
consistent with the virial theorem T+2U = 0, where T, U are respectively kinetic and potential 
energy. This is because the contribution from r ~ R is all potential while r > R is all kinetic. 

Finally with the calculated action we can estimate the lifetime of our universe by using 
eq. fllEp. Using the volume of the past light cone we get consistency if Sb > 500. This is 
easily satisfied, already for X/A ~ 10. 
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Having established the stability of the vacuum on the plateau at the present, one may ask 
the same question for the previous stages of evolution of the universe. We do that by "going 
backward in time" , as later stages of cosmological evolution are obviously more under control. 
So the first question is: what are the conditions for X not to be driven away by thermal 
corrections? Intuitively one might expect that at high temperature T a flat direction X with 
renormalizable interactions is driven to the origin by the positive ~ T 2 corrections to its mass. 
However, this is true only for X < T, as for X 3> T all the modes coupled to X get heavy and 
decouple. In these regions the potential is determined by 1/X suppressed couplings to the 
remaining light fields in the Kahler function. This corrections provide a ~ T 2 /X curvature to 
the flat direction field. This is negligibly small if the maximum temperature of the universe, 
namely the reheating temperature, is much smaller than the X. For a reheat temperature 
^reheat < A 2 X gaugino condensation in the strong group is not affected by the temperature 
corrections and the zero temperature potential is dominant. Note that for values of X to d a y 
favoured by our models the above condition is automatically satisfied when T reheat satisfies the 



gravitino regeneration bound[|4]. We conclude that it is consistent to survive in this local 
minimum during the big bang. The next step behind deals with the evolution of X during 
inflation. Here, however, the theoretical situation becomes less under control, especially in 
the absence of a well defined model of inflation. 



4 Models 

We now illustrate some explicit models built along the ideas of Section 1. We name each of 
them by its Gs x(?b gauge group factor. 

4.1 SU(2) x SU(2) B i x SU(2) B2 - 

The simplest choice for the strong group is SU(2) with a matter content given by four 
doublets (two flavors) Qi,I = 1, ... ,4. In the massless limit this sector has a global 577(4) 
flavor symmetry. Any anomaly free subgroup of it may be weakly gauged and play the role of 
the balancing group Gb- However we find that the simplest way to realize the requirements 
of section 1 is to consider Gb = 577(2) x SU{2) ~ 50(4) C S77(4). For instance, the smaller 
subgroups U(l) or SU(2) would not be asymptotically free. The need for asymptotically free 
balancing group seems to be generic for the mechanism to work. 

The matter content under SU(2) S x SU(2) m x SU(2) B2 x 577(5) is 

Q = (2, 2, 1, 1), Q = (2, 1, 2, 1), R = (1, 2, 2, 1), 
<?3©g2©gi = (1,1,2,3© 2© 1) q 3 ®q 2 ®qi = (1,2, 1,3© 2 ©1) (24) 

where the q^s have been split into their SU(3) x SU(2) x £7(1) irreps. Notice that we add a 
pair of standard model singlets qi, q\ in order to cancel the SU{2) bi x SU{2) B2 global anomaly. 
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The tree superpotential is given by 

W = R (\QQ + A3M3 + A 2 <M 2 + Ai^gi) . (25) 
The classically flat direction of interest is X = (detR)^. Along X, we have 

<R>=X^l °A (26) 

so that SU(2)bi x SU{2) B2 is broken down to the diagonal SU(2) D , and conditions I) and 
II) of Sect. 2 are satisfied. In the regime in which SU{2) D (and the stardard model gauge 
group) is weak, the only relevant non-perturbative dynamics is that of SU{2) S . (We will 
later discusss how good this assumption is.) Then, below the scale of SU(2)s confinement, 
the effective superpotential is W e ff = XXX 2 . The minima of the effective superpotential are 
determined by the RG evolution of the wave function, as discussed in sect. 2. Let us analyze 
this model in more detail. Indicating by fi the renormalization scale we write the RG equations 
(RGE) in terms of r = ^2 ln(/i/M P1 ). Using a vector notation, the gauge (3 functions above 
the scale X are 

d , 1 1 1 1 1 , , 43 N 

tMwa'A'a^^' 1 '- 3 *-^ (27) 

where we have taken g B1 = g B2 = g B , so that the balancing group is 5*0(4). The equations 
for Zx and for the Yukawa couplings are 



dln(Z 



x) 



dr 

din A 2 



dr 



3g% - 2A 2 - 3A^ - 2A^ - \\ (28) 
6A 2 + 3A 2 ; + 2A 2 . + A 2 - 6g% - 3g 2 s (29) 



^1 = 2A 2 + 7A 2 + 2A 2 + A 2 -6^-^ 2 -^ (30) 
= 2A 2 + 3A 2 + 6A 2 + A 2 -6^-3^ 2 -^ 2 (31) 

din A 2 



dr 



2A 2 + 3A^ + 2A 2 : + 5A 2 - 6g 2 B (32) 

(33) 



where Aj's are the running couplings (e.g. A 2 = X 2 /(Z x ZqZq)). By identifying \x with X, the 
minima of the potential are given by the solutions to dZx/dr = with d 2 Zx/dr 2 < 0. The 
quantity Z x is readily computed from eq. (|28|), with the help of the other RGE. Notice to this 
effect that the only terms that can make Zx > are the gauge contributions to the Yukawa 
RGE. These contributions are indeed big for A 2 , and they are sizeable also for A 2 and A|. 
Therefore for small A 2 , the condition Zx < cannot be satisfied. On the other hand, when A 2 
is bigger than the other Yukawa couplings, one gets Z x < 0. Still, in this regime, at a given X, 
A 2 cannot be made too big, otherwise it hits a Landau pole below M Pi . Such considerations 
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Figure 2: The region of parameters where the SU(2) x SU(2) x SU(2) model can have a 
local supersymmetry-breaking vacuum at (X) = 10 10 GeV. 



determine the regions of parameter space where our mechanism works. In fig. 1 we display, 
for different values of A, the region in the plane (<7b,Ai) for which there is a minimum at 
X = 10 10 GeV. We have assumed A 3 = A 2 to reduce the number of parameters and because 
A3 ~ A2 is motivated by GUT considerations. The boundaries of the region in the figure are 
explained as follows: i) the right one is given by Zx < 0; ii) the top one is essentially dictated 
by the absence of Landau poles in Af below the Planck scale; iii) the left one is determined 
by A^ < 3g% from Z x = 0, see eq. (ffip. So it is clear that there is a significant region of 
parameter space where this model behaves, in the observable MSSM sector, as a conventional 
gauge mediated supersymmetry breaking model, with messenger scale around 10 10 GeV. We 
remark that this model shares similarity with an interesting proposal by Izawa and Yanagida 
23] and by Intriligator and Thomas [p4 |, where the role of our R is essentially played by a 
set of gauge singlets. The phenomenological version of that model, mentioned in Ref. |24j 
and discussed at length in Ref. @], is unfortunately not calculable. In that model, like in 
ours, one must live in a local minimum, as the QQ and the qq restore supersymmetry at 
X = 0. However this local minimum, if it exists, is presumambly at X < A, as the Goldstino 
superfield X is a gauge singlet and the perturbative superpotential drives it towards the 
origin (the perturbative corrections to K(XX^) were ignored in Ref. |24|, |3||). In this region 
the Kahler potential is not calculable and one has just to assume: 1) the occurrence of such a 
local minimum; 2) its quantum stability, as there is no obvious parametric suppression of the 
tunnelling rate at X ~ A. 
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Around the minimum the real part of X gets a mass 

which is parametrically similar to the sfermion masses in eq. ([13]). At this stage, the imaginary 
part of X is the R-axion. Notice that X, the axion scale, as far as gauge mediation is concerned, 
can consistently be in the 10 9 — 10 12 GeV axion window P3[ . So X also provides and interesting 
QCD axion candidate. However there is still the possibility that ImX gets a mass from direct 
sources of explicit R breaking. One of these can arise from the mechanism that cancels the 
cosmological constant f2~6|| . Another is the SU(2)p, anomaly, but this is negligible when the 
balancing group is weak (see the comments below). 

Actually also SU (2) D has a low energy dynamics that leads to an additional term Wp> ~ 
A#X 2 , where Ap = A12 is the microscopic strong scale of the balancing 577(2) factors. 
The combination of the two gaugino condensates restores supersymmetry, as the low energy 
W ~ A 2 X + AbX 2 has no effective R-symmetry. However, since we are interested in local 
minima we are not worried by that, as long as Ap is small enough that the supersymmetry 
preserving minimum is very distant. This is indeed what happens for the values of gp{X) 
that are allowed in Fig. 1. Notice that the relative shift in the vev is SX/X ~ 87T 2 Ap/mx, 
which for gp(X) < 1 is less than 10 -20 , while another supersymmetry preserving minimum 
(in addition to those at X = 0) exists only at X > Mpj! 

Some more comments are in order. These concern the range of X where the mecha- 
nism works. In particular X cannot be too close to A, otherwise our perturbative approx- 
imation to the Kahler potential breaks down, due to non-perturbative corrections from the 
SU (2)s dynamics. Indeed at A > A we expect the non-perturbative corrections to K(X, X*), 
to be dominated by the lowest term in the operator product expansion. This is given by 
/ci 4 #(iy 2 )(W/ 2 )/(XXt) 2 , whose coefficient is ~ ^2 as it is determined by a 1-loop diagram 
involving the heavy Q,Q and four external vector superfields. After substituting the gaugino 
condensate the above operator leads to a correction SK(X, X') ~ 87r 2 A 4 /(A 2 XX^). Thus the 
effects induced on the Kahler curvature vanish like (A/X) 4 . Using the above formula for m 2 x , 
we estimate that this non-perturbative effect can be reasonably neglected for X > 10 9 GeV. 

At the other end of the range, for large X, operators induced by Planck scale phisics may 
affect the flatness of the superpotential. We find that, barring the mass term W = MpiX 2 , 
the lowest correction X 4 /Mpi is tolerable when the minimum is at X < 10 11 GeV. For higher 
messenger scales this term will also have to be suppressed. We notice that while the X 2 and X 4 
operators are in principle allowed by the symmetries of the model, it is not totally unreasonable 
to assume they are absent. This is because in supersymmetric theories the absence of a given 
term in an effective superpotential (now eq. fl2"o] ) itself) does not necessarily have a symmetry 
explanation in the low-energy theory. (This is for instance seen in examples of strong gauge 
dynamics.) 

Before concluding we stress that the phenomenology of this model is that of conventional 
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gauge mediation with two families of messengers in the fundamental of S77(5). One could 
consider generalizations based on Sp(m) x Sp(m) x SU(2m), where m = 1 is just the original 
model. These work essentially in the same way, but have 2m families of messengers. 



4.2 SU(2n) x SU(2n) with adjoint matter. 

Another possibility is to consider Gs = SU(N) S and Gb = SU(N) B , with even N = 2n 
and with the following matter content under SU(N)s x SU{2)b x S77(5) 



Q = (N,N,1), Q = (N,N,l), #=(1, AT 2 -1,1), g = (1, 7\T, 5), g = (l,A>,5). (35) 



where the entries represent n x n blocks. To motivate this choice of model and W we remind 
that in the absence of a superpotential, there exists an iV — 1 dimensional family of flat 
directions parametrized by iV — 1 invariants built from R: X p = (TiR p ) 1 ^ p with p = 2...N. 
At a generic point on this vacuum manifold, SU(N)b is broken down to [/(l)^ 1 and there 
are just N — 1 massless chiral and N — 1 massless gauge superfields. There are, however, 
enhanced symmetry directions (lines) along which a non-Abelian subgroup is unbroken. Here 
the number of massless gauge and chiral superfields is bigger. A suitable superpotential can 
lift all but one flat direction. For SU(mn) groups (to, n integers) adding W = TrR m+1 , lifts all 
moduli but u m = TiR m , along which SU(mn) — > SU(n) m . The special thing about SU(2n) 
groups is that this job is done by a trilinear invariant Tri? 3 . This motivates our choice of even 
N with W given by eq. fl36|) . 

Now, along X, see eq. (j37D, the balancing group is broken down to SU{n) x SU{n). All 
other components of R as well as Q, Q and q, q have a mass of order X and decouple. The 
low energy dynamics of SU(N)s generates a linear superpotential W = XX A 2 . One can go 
through the same analysis we applied to SU(2) 3 . In the models with iV = 2n < 4 the balancing 
group is not asymptotically free, and this makes it impossible to satisfy the minimum condition 
Zx < 0. Moreover models with iV > 10 have too many messengers, so that Landau poles in 
the evolution of the SM gauge coupling are hit before the GUT scale, unless X is very close to 
Mq, bringing partially back the supersymmetric flavor problem. We find that a model which 
works very well over a wide range of scales is 577(6) x SU(6). Let us discuss it in more detail. 
The gauge coupling RGE are 



The classical superpotential is 




(36) 



Again, at the classical level, there is a flat direction along which only R ^ 




(37) 



d , 1 1 1 1 1 



) = (1,12, -3, -7, 
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(38) 
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Figure 3: The region of parameters where the SU(6) x SU(6) model can have a local 
supersymmetry-breaking vacuum at (X) = 10 10 GeV. 



The equations for Zx and for the Yukawa couplings are 

(39) 

(40) 

(41) 

(42) 

(43) 
(44) 

where again the barred Yukawas indicate running couplings. 

In Fig. 2 we show, for different values of A the region in the (gB,k) plane where the model 
leads to a stable minimum at X = 10 10 GeV. Again we assumed A2 = A3. The two disjoint 
regions arise as follows. The thin region around the line gs = (2/3) k, corresponds to the 
points where Zx = is saturated essentially by k, see eq. fl3"9"|). To the immediate right of 
the thin area, A2,3 are small but non zero. This is enough for the negative contribution of g% 
and #2 to their RGE to make Zx positive. This excludes this intermediate area. However 
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when A2,3 become large enough, the positive terms in their RGE dominate and Zx becomes 
negative again. This explains the big allowed band to the right. However as ^2,3 are further 
increased, they hit a Landau pole below Mpi, and the region of large qb is again excluded. 
The range where the model works becomes even wider as the scale X is increased above 10 10 
GeV, as the constraint from Landau poles gets weaker. Also, as for the previous model, the 
effects on W of the low energy SU(3) 2 are negligible as soon as #b(A) < 1. 

4.3 SU(5) 3 . A special model. 

As a special case of our general idea one can consider the possibility of identifying Q, Q 
with q, q. In other words the messengers themselves trigger supersymmetry breaking. This 
can be achieved in a model with full gauge group 577(5) i x SU(5) 2 X SU(5) 3 and with matter 
content given by R(l, 5, 5), Q{5, 1, 5) and Q(5, 5, 1). For phenomenological applications we 
are interested in a variant (low-energy limit) of this model, one in which S77(5)3 is already 
broken to SU(3) x SU(2) x U(l) and the standard matter transforms in 5© 10 of SU(5)s. To 
simplify the discussion let us first focus on S77(5) 3 . In view of its role we treat SU(5)3 as a 
weakly gauged factor, and neglect possible effects coming from its non-perturbative dynamics. 

The space of D-flat directions of the theory is parametrized by 3 baryonic invariants 
u = det(i?), v = det(Q), z = det(Q) plus a meson Y = RQQ. A renormalizable superpotential 
consists just of the term 

W = XRQQ = XY. (45) 

This superpotential lifts most flat directions. The classical equations of motion, 8rW = 0, 
8qW = and dqW = 0, when contracted to form invariants give the constraints 

uv = uz = vz = Y = 0. (46) 

Thus the space of classically flat directions consists of three one- dimensional baryonic branches 
u,v,z with one common point at the origin. Notice that this result holds also for the model 
in which one gauge factor, say SU(5) 3 , is not gaugedF]. 

The class of models to which this one belongs have been discussed in ref. |27J where it 
was shown that in the presence of the trilinear superpotential (and treating SU(5)3 as weakly 
gauged) they break supersymmetry at any finite point in field space. Furthermore there is 
a sloping direction z — > oo along which the vacuum energy becomes asymptotically zero. 
However, the dynamics which interests us takes places far away along the u (or v for that 
matter, since SU(5)i and SU(5) 2 can be interchanged) flat direction. It is just a special case 
of the dynamics discussed in Section 2. Along u, R = X X 1, so that u = X 5 . This breaks 
the factor SU (5)2 x S77(5)3 down to the diagonal SU(5)d, which will be later identified with 

2 In that case the space of D-fiat directions consists of the same baryons plus a meson matrix = (RQQ) 1 ^, 
where A, B are SU (5)3 indices, supplemented by the classical constraint detF = uvz. Again the superpotential 
is W = ATrF, and the condition of F-flatness leads to a vanishing meson matrix = Qb9q a W — XYg and 
hence back to eq. (f46[). Thus even for a global SU(5)s only the baryonic branches of moduli space survive. 
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the standard model gauge interactions. Apart from X, all the fields in R, Q, Q get mass. 
Thus the SU{5) i gaugino condensation generates a superpotential W e ff = A| = AXA^. A 
few comments are in order at this point. First of all, this is now an exact result. There are 
no corrections to it coming from the SU (5)2 dynamics. This is not totally surprising, as that 
group is Higgsed. This result can also be established by studying the global symmetries of the 
microscopic Lagrangian and using holomorphy. In particular the baryon number U(1)q with 
charges (0, 0, 1, —1, 0, 5) for (A, X, Q, Q, A\°, A^ ), rules out any dependence of the low energy 
superpotential on A 2 . 

Now, let us analyze quantitatively the realistic scenario where SU (5)3 — > SU (3)' x SU(2)' x 
U(l)'. To avoid confusion with gauge group suffixes we rename SU (5)2 — > SU (5) b (balancing) 
and SU(5)i — > SU(5) S (strong). The previous discussion is easily generalized. Now R — > 
R2 © R3 = (1, 5, 2 © 3) and similarly Q — > Q 2 + Q3. The flat direction analysis is unchanged. 
The vev of R breaks SU(5) B x SU(3)' x SU{2)' x U(l)' -> SU(3) x SU(2) x U(l), so that 
the low energy gauge couplings are 

1 1 1 

+ 327^ (47) 



gKx) g 2 B (x) g ?{x) 

where the notation is obvious. Notice that at the scale X the gauge couplings gi are modified 
in an SU (5) invariant manner, so that gauge unification is safe. However, the coupling above 
X will be somewhat bigger (like y/2). This we find to be the major source of constraints for 
this model, as we discuss in more detail below. 

The SU(5)\ nonperturbative dynamics goes along unmodified. The field R has however 
broken up into two pieces i? 2 ,3 whose wave function terms evolve differently and which we 
indicate by Zx 2 and Zx 3 - Now, along the flat direction R 3 = X 3 l3 X 3, R 2 = A 2 l 2 x2 with 



X3/X2 = y Zx 2 /Zx 3 in order to have vanishing SU(5)b D-terms. As before, we parametrize 
the flat direction with X 5 = A|X| = deti?. The Kahler potential along X is now given by 

K(X,X*) = hZl 2 Zl 3 XX* (48) 
so that the minimum condition is 

5 |£= 2 |^ + 3§£i = 0. (49) 

6x 2 ^X 3 

The gauge coupling RGE are 



d , 1 1 1 1 1 x „ 58, 

Tr { Aw s 'Trww ) = {w ' w '- 2 '- 6 '~ ] - (50) 

while for the Yukawa couplings we have 



d^ n ( Z x) 24 2 3 /2 8 , 2 1 /2 ^ 2 o\2 /r-,\ 
^ = —9 B + -^92 + -^93 + q01 - 2A 2 - 3A 3 ( 51 ) 
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^1 = 13 A| + 2AS-f 9 l-f 9 l-f 9? -l 9f (52) 

^1 = 12A1 + 2AI - f ^ - f ff I - 3 ff ? - (53) 

(54) 

This model has no problems with Zx becoming positive. This a consequence of 577(5) b being 
"strongly" asymptotically free. In a situation like this, the minima in the 1/Zx potential are 
well pronounced. In fact they are so pronounced that the parameter space of this model is 
strongly constrained by the avoidance of Landau poles in A|, A|. The constraint is made even 
more acute by eq. fl4T|), which implies that both g B and gf have to be larger, in practice by 
about a factor of 2, than the effective gf. This implies that the Yukawas that saturate Zx = 
have to be rather large too, and their evolution above X quickly reaches a Landau pole. We 
find that the request of no Landau pole below 10 17 GeV and of perturbative gauge unification 
sets a lower bound on the scale X 

X > 10 14 GeV (55) 

For such an X we find solutions with the following values of the parameters: g%{X) ~ 1 — 1.2 
and A|(M(?) — ^K-^g) ~ 1.8 — 3 and gf^jT < 2. This is still in the perturbative domain. 
However, by lowering X we immediately get large Yukawas at the GUT scale. Notice that for 
X ~ 10 15 GeV, we start recovering the flavour problem because of gravity contributions to soft 
masses, thus weakening the motivation for gauge mediated supersymmetry breaking. So we 
can say that this model works in a fairly limited but still non vanishing range X ~ 10 14 — 10 15 
GeV. 

Again we should comment on all additional contributions that may affect the scalar poten- 
tial. In this model, given the small ratio A/X ~ 10 -5 , the non-perturbative corrections to the 
Kahler metric from the SU(5)s dynamics are clearly negligible. But we are very close to the 
Planck scale. The lowest operator that can affect the flatness is SW = detX/M Pl ~ X 5 /M Pl . 
Then it is clear that for X ~ 10 14 GeV this correction must be absent, or be very suppressed. 
If this is the case, flatness is preserved, as the higher order terms like detX 2 / M 7 Pl are very 
suppressed. One last source of perturbation can arise from whatever mechanism relaxes to 
zero the cosmological constant by providing an effectively constant term ~ A 2 Mpi in W. This 
term creates in lowest order a linear potential 5V ~ A 4 ReA jMp\. This leads to a shift SX in 
the vev of X, which is 

where the (8vT 2 ) 2 arises from the expression for m x , see eq. ([34]). Then again, X ~ 10 14 GeV 
is still acceptable, while X ~ 10 15 GeV starts being problematic. This bound is amusingly 
similar to the one provided by the flavor problem. 

In spite of its very constrained range of validity we think this model is an attractive example 
of a very elegant way of generating and communicating supersymmetry breaking. Are there 
any special phenomenological signature of models like this? The answer is in principle yes. 
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Indeed in this theory the messengers of supersymmetry breaking are not exhausted by the 
chiral superfields Q,Q. There are also the massive vector superfields in the coset SU(5)b x 
SU(3)' x SU(2)' x U(1)'/SU(3) x SU(2) x U(l). These get a mass from X and feel the effect 
of Fx 7^ 0. A set of these massive vectors transforms as the adjoint of SU(3) x SU(2) x U(l) 
and couples directly to standard matter. We may call these heavy objects the axial gluons, 
W's and B. Indeed already at one loop they affect the sparticle parameters. Since the field 
X is the only source of infrared cut-off, the form of the correction is fixed by the 1-loop wave 
function renormalization, as for K(X,X') is Section 2 



where the notation is the same as in eq. (|T|). Notice the interesting fact that this 1-loop cor- 
rection generates only A terms at order a/ Arc. The lowest contribution to the chiral preserving 
soft masses is zero as dxd x \ ln(JX*) = 0. Chiral preserving mass squared are generated at 
0(a 2 ) by both the above 1-loop operator and by a host of new 2-loop effects. In this way the 
correct scaling of gaugino and sfermion masses is preserved. This cancellation of dangerous 
1-loop effects is the same as the one observed in ref. p8| , |29| for the case of standard matter 
directly coupled to the messengers via Yukawa couplings. We stress once more that in the 
superfield language the cancellation simply follows from the coincidence of the infrared regu- 
lator X and the source of supersymmetry breaking. It does not work, for instance, when the 
particles in the loop get mass from two independent spurion superfields X and Y for which 
X/Y is not a c- number. As it was shown in ref. [^, the cancellation works only for 
the leading contribution in an expansion in powers of (A/X) 2 . The sub leading terms corre- 
spond to higher derivative terms in the effective action for X and SM matter, which are not 
determined by the RG log and for which the above argument does not apply. These higher 
order effects are however unimportant for this model as A/X <C 1. A more detailed analysis 
is needed to determine the effects on the sparticle spectroscopy of the terms induced by the 
axial vector superfields. It would be interesting to find experimentally important differences. 
This we leave for further study. One sure thing is that A terms are now induced directly at the 
messenger scale and they are of the same order as gaugino masses . This may have important 
implications for the case in which the fi term in the Higgs sector is provided by a singlet S 
with couplings W = S 3 + SHiH2- In conventional gauge mediated models the smallness of the 
RG induced A terms in this sector makes the correct breaking of the electroweak group more 
difficult. This is the case even when the messenger scale is fairly high. It would be interesting 
to study the impact of the new boundary conditions. 

5 Upper Limits on the Messenger Scale from Nucle- 
osynthesis 

The messenger scale is M is constrained by laboratory and astrophysical considerations. The 
first is the supersymmetric flavor problem, which is a primary motivation for gauge mediated 




(57) 
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theories. It costrains the messenger scale to be less than about 10 15 GeV 0. 

The purpose of this section is to introduce a stronger constraint coming from the require- 
ment that the NLSP decays do not spoil the successful predictions of primordial nucleosyn- 
thesis. These constraints will turn out to be fairly independent of the nature of the NLSP 
which is either a neutralino or the right handed stau depending on the number of messengers 
iV as well as the messenger mass M as shown in fig. 3. Let us first discuss the case when the 
NLSP is a neutralino, x°- The decays of x° are studied in detail in reference []3~0[| . There it is 
shown that, in general, the dominant decay modes of x° are into a photon and gravitino and, 
if kinematically allowed, to a lesser extent into Z° and gravitino. The decay rates are given 
by: 

5 

r( x °^ 7 G)^cos 2 ^ T ^ (58) 

r( X ° - ZG) * sm 2 6 w -^- 2 (l - (59) 

16ttF 2 \ m x° / 

The photon of the dominant decay mode can affect nucleosynthesis because it can break- 
up heavier nuclei into lighter ones [32], [33]]. The "damage" that these decays can do 
to nucleosynthesis depends on the product of the mass of the decaying particle times its 
abundance relative to baryons (or photons): 

Tl 

d 1 = m x o — . (60) 

This quantity is somewhat model dependent, as it depends on the sparticle spectrum. For 
example, the lightest neutralino in typical gauge mediated theories is mostly a B-ino and its 
abundance relative to baryons depends on its self-annihilation rate and consequently on the 
slepton spectrum. This in turn depends on N and M and other features of the messenger 
sector [fjEH - Nevertheless, the order of magnitude of d 1 is easy to estimate provided that the 
typical mass scale of the light supersymmetric particles is about 100 GeV, resulting in a self- 
annihilation cross section of order 2w(a/m)' 2 ~ 10~ 35 cm 2 mj~ 2 . Here mioo is the typical mass 
of light sparticles in units of 100 GeV. This yields a mass times abundance per unit entropy 
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= m— ~3x 10~ n GeV, (61) 
s 



corresponding to g? 7 ~ 0.3 m 10 o GeV. This implies that [[32], |33| as long as the lifetime of 
the decaying particle is less than 10 7 sec the damage that the photonic decays can do to 
nucleosythesis is negligible. They imply an upper limit of the scale of supersymmetry breaking 
\J~F < 10 10 GeV or equivalently 

M < Nm% 2 10 15 GeV, (62) 

where N is the number of messengers. Therefore we see that the constraint from photode- 
struction of light elements is quantitatively similar to that from the flavor problem. 
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Figure 4: The regions in the N-M parameter space where the neutralino or the right- 
handed stau (f) are the next-to-lightest supersymmetric particle (NLSP). Above the upper 
curve the gauge coupling constants become non-perturbative before the GUT scale. 



A stronger constraint arises from the hadronic decays of x° @ • These occassionally lead 
to energetic nucleons which subsequently thermalize by colliding with the ambient protons 
and 4 He and in the process they create a "hot line" along which a one dimensional version of 
nucleosynthesis can take place, potentially resulting in overproduction of the light nuclei D, 
3 He, 6 Li, 7 Li [R2 . D and 3 He arise mostly from "hadrodestruction" - an energetic nucleon 



breaks up an ambient 4 He nucleus into 3 He or D. 7 Li and 6 Li arise from "hadrosynthesis" 
in which 3 He or T or 4 He are sufficiently activated along the "hot line" to overcome their 
Coulomb barrier and synthesize 6 Li or 7 Li nuclei. The analysis is fairly involved and the 
reader is referred to reference p3J for the details. The crucial quantity measuring the 



magnitude of the "damage" of hadronic decays on the primordial nuclear abundances is 

d B = d 7 r* B , (63) 

where r* B is essentially the hadronic branching ratio r B of the decaying particle - it is defined 
carefully in references p4| , 32| . The photon in the decay x° can convert to quarks 



with probability of order a, and therefore r B can never be smaller than about 10 2 . The 
contribution to the hadronic branching ratio from the decay x° ~^ ZG is 



r B ~ 0.7 tan 2 ^ 1 f) , (64) 

where 0.7 is the hadronic branching ratio of the Z. This is typically smaller than 10~ 2 as 
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long as the x° mass is less than about 150 GeV because of phase space suppression. So from 
now on we take r# to be its minimal value 10~ 2 ; if it is larger then our arguments will get 
stronger. The value of the corresponding ds is 3 x 10 -3 GeV. From figures 5 of reference 
|32|j we see that even for such a small value of ds the hadronic showers lead to an extreme 



overproduction of 7 Li, 6 Li and D + 3 He, if the lifetime of the decaying particle is longer than 
10 4 sec. This shows that the lifetime of x° must be significantly smaller than 10 4 sec. How 
short the lifetime has to be before we are safe is not yet known to us. For lifetimes less than 
10 4 sec we are not aware of a complete analysis - including the crucial hadronic jets - of 
the effects of late decays on all the light element abundances, in particular 6 Li and 7 Li. One 
complication is that ordinary nucleosynthesis is occuring simultaneously with the decays. For 
this reasons it is not possible yet to quote a precise upper limit on the lifetime of a decaying 
particle from nucleosynthesis. A safe upper limit for the lifetime is about 1 sec, which marks 
the beginning of nucleosynthesis. This leads to an upper limit of about 10 8 GeV for the scale 
of supersymmetry breaking \/F or, 

M < Nm%l 10 11 GeV. (65) 

This is a safe upper limit on the messenger scale and is significantly stronger than the limit 
of M < 10 15 GeV coming from the flavor problem. However, it may be too strict and a 
detailed analysis may allow for somewhat longer lifetimes [571]. Thankfully, the upper limit on 



M scales as the square root of the lifetime, so that even if lifetimes of order 100 seconds were 
allowed they would only increase M by an order of magnitude. This would still lead to an 
upper limit on M much stronger than the limit coming from the flavor problem. Since most 
nuclear abundances, in particular that of 4 He, are in full bloom by the end of 180 seconds we 
expect that lifetimes longer than 180 seconds are excluded. The point is that the damaging 
processes of hadrodestruction and hadropro duct ion - which lead to the overabundances of 
D+ 3 He and Lithium respectively - rely mainly on the existence of ambient 4 He which is 
completely formed by 180 seconds. Furthermore, the universe is sufficiently diluted by 180 
seconds that a reduction of the excess abundance is unlikely. 

When the NLSP is a right-handed slepton the upper limit to the messenger scale should 
be even smaller than in eq. (|65|). The point is that the lightest right handed slepton is the 
right handed stau which decays into a tau and a Goldstino. Since the tau has a large hadronic 
branching ratio, the hadronic damage factor ds becomes 0.2 GeV which significantly exceeds 
that of x°- 



The upper limit in eq. (|65"D has significant implications on model building. It disfavors 
theories in which the messenger scale results from the competition between renormalizable and 
non-renormalizable Planck-supressed operators, since these result in large messenger scales. It 
favors theories in which the messenger scale is the result of renormalizable competing forces, 
as is the case here. Even then, it excludes our most favored theory, namely SU(5) 3 . 



One, perhaps unappealing, way out of the constraint in eq. (p5|) is to postulate some 
degree of i?-breaking which can result in a microscopic lifetime for the NLSP which is totally 
decoupled from supersymmetry breaking. 
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Finally, notice that the longest NLSP lifetime consistent with the flavor problem (M < 10 15 
GeV) is 10 7 sec. This implies that the NLSP decays cannot have an effect on the spectrum of 
the microwave and diffuse background radiation |33|. [36| . 

It is interesting that the nucleosythesis constraint implies that the scale of supersymmetry 
breaking breaking has to be in the "narrow" range between 10 4 and 10 8 GeV. To avoid the 
gravitino problem in the high end of this range - between 10 6 and 10 8 GeV - it is necessary 
that the reheating temperature of the universe be relatively low. 



6 Conclusions 

Gauge mediated supersymmetry breaking provides a very attractive alternative to gravity 
mediated scenarios. In gauge mediation the only source of flavor violation which is relevant in 
low-energy phenomenology is given by the Yukawa matrices of standard quarks and leptons. 
This leads to the supersymmetric generalization of the GIM mechanism and to the distinctive 
experimental signature that no major deviation from the SM is expected in flavor and CP 
violating processes. This is consistent with our present experimental information, and is the 
main motivation for gauge mediation. The gravity scenario, on the other hand, has a generic 
difficulty with flavor violation. Nonetheless, the models that can be made consistent with the 
present bounds (eg. by using horizontal symmetries) still predict that some major deviation 
from the SM (like \x — > ej or an electric dipole moment of an elementary particle {40|) will 
sooner or later be discovered. Whether gravity or gauge interactions mediate soft terms may 
be answered by accurate low energy experiments. 

The breaking of supersymmetry and its mediation (via gauge interactions) at scales much 
below Mpi should be well described by field theory, with no active role played by quantum 
gravity. This motivates the study of theories where supersymmetry is broken by some non- 
perturbative gauge dynamics. This allows for an attractive field theoretic explanation of the 
smallness of the weak scale itself, not just its quantum stability. It is clear that the discovery of 
compelling field theory models with dynamically driven and gauge mediated supersymmetry 
breaking would give this scenario another non-trivial boost over the gravity mediated case. 

These facts are at the basis of numerous recent attempts to build realistic theories wedding 
dynamical SUSY breaking to gauge mediation. The main point of our investigation has been 
that many of the aesthetic and phenomenological disadvantages of previous attempts are easily 
avoided by relaxing the request that supersymmetry be broken in the true vacuum. We have 
characterized a class of theories for which the vacuum where we live is typically a false one, 
located on a supersymmetry breaking potential energy plateau, far away from from other 
energetically more favourable minima. This last property can be achieved quite naturally by 
means of Witten's inverse hierarchy phenomenon and endows the false vacuum with a lifetime 
which is much longer than the age of the Universe. This makes it in practice a stable vacuum. 

We have characterized a class of calculable theories with a fairly simple structure. The 
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breaking of supersymmetry is communicated to the MSSM sparticles in a minimal number 
of steps. A measure of this simplicity is the fact that the Goldstino superfield itself gives a 
mass to the messengers. This is not the case, for instance, in the standard scenario by Dine 
et al. |2j, which contains an additional sector and where the Goldstino is not directly coupled 
to the messengers. Those models have no supersymmetric ground state, but this does not 
lead to any clear advantage as the global minimum breaks color and SU(2)w at a high scale. 
Thus the messenger sector has to live in a local minimum. It seems that in order to have 
a phenomenologically acceptable true minimum considerable further complication is needed 
|4j]. Notice that these minima are the analog of those that exist in our models, see eq. (1), 
at X = 0. Indeed in gauge mediation with a large messenger mass there can yet be other 



deeper minima along the flat directions H2 = Li [^8] of the MSSM. This is because in gauge 
mediation the soft mass m^ 2 + m| along this flat direction goes negative very quickly due 
to the large stop contribution [f£|. So it seems that in one way or another gauge mediation 
brings us the idea that we live in a false vacuum. 

Attempts to build theories with just two sectors have flourished in the last year. In 
particular the models proposed in Ref. ,0 have a close resemblance to our SU(5) 3 example. 
Those models are however phenomenologically non viable, as squarks masses turn out to be 
negative. This is due to the presence of relatively "light" matter supefields with positive 
and large supertrace. Since the potential is stabilized by 1/Mpj operators, the scale where 
soft terms (and the supertrace) are generated is fairly large ( > 10 15 GeV, barely compatible 
with the flavor problem), the positive suprtrace leads to large and negative squark masses via 
two loop RG mixing. We stress that our SU(5) 3 model bypasses the problem by having no 
such light matter: everything gets a ~ X mass via just one Yukawa coupling. Other new 
mechanisms were explored in Refs. |7], |8| , where the role of X is played by a composite gauge 
invariant field. Those models, however, involve massive parameters in the superpotential, so 
they do not yet represent a full microscopic theory. 

In previous cases the messenger scale was either fairly low ~ 10 5 GeV, for renormalizable 
theories, or fairly high > 10 15 , for non-renormalizable ones. In our scenario we have a double 
dimensional transmutation, one for Fx (gaugino condensation) and one for X (Coleman- 
Weinberg mechanism). By that we can cover the continuum of theories where messengers have 
a mass that varies from 10 8 ~ 9 to 10 14 GeV, and where 10 6 ^ 7 GeV < F x < 10 9 GeV. Models 
that live in this intermediate region can more easily avoid the gravitino problem of usual 



renormalizable models and gravity has no observable flavor effect. Indeed we also pointed 
out that models with a large messenger scale > 10 12 GeV have serious problems with standard 
Big Bang nucleosynthesis. This still keeps a large slice of the continuum unconstrained. 
Models with X > 10 12 GeV, including the otherwise very compelling S77(5) 3 , are in need of 
additional interactions which allow the NLSP to decay in less than a second. It seems that a 
possibility would be to add a small source of i?-parity violation, but with the aesthetic and 
logical (remember flavor physics) problem that this addition entails. Another totally different 



way out is to assume a period of GeV scale inflation 42], with GeV scale baryogenesis. 



A possible direction for the future is to further enlarge the taxonomy of theories with 
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gauge mediated supersymmetry breaking. For instance, one hope is that with a large pool to 
choose from, a still unresolved problem like the origin of fi may be solved in a natural way |[28|| . 
Indeed the boundary conditions on soft terms plays a relevant role in this issue, making, for 
instance the solution by the addition of a singlet S more or less attractive. A crucial aspect 
is the size of A terms. We have seen that models with heavy vector supefield messengers can 
lead to a rather different boundary condition with respect to the minimal case, where A terms 
are generated by gaugino masses via RG evolution. 
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Note Added 

As this paper was about to be submitted a very interesting paper by H. Murayama, |hep- 



ph/9705271| , appeared. In that paper a model similar to our SU (5) 3 is discussed. 
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Appendix 



In this brief Appendix we discuss the dynamics along the flat direction X as the field VEV 
becames of order A. As an example we choose to illustrate the SU(5) 3 model, because it is a 
very attractive one. Nonetheless, in the other cases, the essence of the discussion is the same. 

As already said the strong dynamics of SU(3)' x SU(2)' x U(l)' can be neglected. The 
full effective superpotential of the confined S77(5)i x SU(5) 2 theory is given by f2"7|| 

W conf = XTrY + A (det(Y) - uvz + A\°u - A£°v) (66) 



where A is a Lagrange multiplier enforcing a quantum modified constraint ]13|. At least for 



A <C 1 we expect the above fields and W con f to describe the configurations of lowest energy. 
Notice that the above superpotential admits no supersymmetric vacua at finite field values. 
Indeed Fy = requires A, Y ^ after which F z = implies uv = 0. This last requirement 
is inconsistent with F u = F v = 0, unless z — > oo. Far away along z all other baryons and 
meson are massive, as it is already shown by the classical analysis, and integrating them out 
yields the effective potential W e ff = (A^A 2 ,) / z 1 ^ 5 f27|. So supersymmetry is asymptotically 



recovered only as z — > oo: this theory has no global vacuum. We are however interested in the 
u branch and in what happens as X is lowered to be ~ Ai. In order to so, we first integrate 
out the meson Y and the field A. This is sensible as they are massive over the entire moduli 
space. As a matter of fact the baryons u, v, z are the only fields that can become very light 
on some branch of the moduli space. We thus obtain 



i 



W ef f = A {uvz + uA\° - vAl°) 5 (67) 

The above equations contains all the interesting information regarding the dynamics as we 
move along X = (w) 1//5 . For a ^ we can make the reparametrization z = z — A^/u. Then 
performing the suitable rescaling of the baryons v , z their mass matrix is roughly given by 



£"mass — (v*, Z*) I ,. 2 \2 \ y\2 )( ~ )• (68) 




where the diagonal terms come from \F V ^ | 2 , while the off diagonal from ji^xl 2 - At large X the 
baryons are massive with zero vev, so that by decoupling them we get back the linear W e ff. 
However, when X is decreased below some critical value ~ Ai the baryons become tachionic 
and develop a vev. This indicates that for X ~ A ly /A we enter the basin of attraction of the 
z flat direction and then slide to infinity. The only point that is crucial for our discussion, 
however, is that no instability develops for X > Ax/A, and thus to tunnel to the asymptotically 
supersymmetric vacuum the bounce configuration has to interpolate between field values of 
order X to values of order Ai/A. 
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